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摘要 :基于 1986 年 到 2011 年 的 Landsat 影像 


统 (LEDAPS) 预 处 理 生成 标准 的 地 面 反 出 


E SUELE: 


,以 南方 人 工 林 分 布 区 域 广 东 省 佛 冈 县 为 例 ,运用 Landsat 生态 系统 自 适 应 处 理 系 
E Landsat 时 间 序 列 堆栈 (LTSS) 用 于 LandTrendr 算法 监测 人 工 林 森林 干扰 


与 恢复 的 长 时 间 序 列 变化 ,分 析 了 连续 24a 森林 干扰 的 年 份 变化 .干扰 量 以 及 干扰 持续 的 时 间 ,验证 了 算法 识别 干扰 的 精度 ,并 


探讨 了 人 工 林 干扰 的 驱动 力 。 结 果 表 明 人 
和 2009 年 的 干扰 面积 均 超过 2000 hm? ,其 中 1987 ,2007 年 两 年 的 干扰 面积 达到 6000 hm? 以 上 。 相 比 森 林 干 扰 的 变化 ， 


四 县 的 森林 干扰 较为 剧烈 ,一般 都 在 1000 hm? , 。 而 1987 ,2002 .2004 .2005 、2006 


\、2007 


四 县 


的 森林 恢复 面积 随时 间 的 变化 相对 平稳 。 通 过 对 佛 内 县 森林 干扰 和 恢复 面积 的 趋势 分 析 , 发现 20 世纪 80 年 代 末 到 90 ERR 


林 干 扰 和 恢复 的 函 
渐 上 升 ,总 体面 积 变化 趋势 高 了 
约 38% ,持续 2a 


的 森林 干扰 和 恢复 的 面积 在 人 
后 ,急剧 干扰 的 面积 远大 于 持续 干扰 ,最 高 约 达 2800 hm ,但 


积 基 本 少 于 2000 年 以 后 的 变化 面积 ,变化 趋势 比 2000 年 以 后 的 显得 平缓 ;从 2000 年 开始 ,森林 干扰 多 
F 森 林 的 恢复 ,但 森林 的 恢复 面积 仍 有 所 提升 。 其 中 , 佛 冈 县 的 森林 干扰 持续 1a 时 间 的 面积 比例 
时 间 约 28% ,持续 3a 时 间 约 25% ,持续 4a 时 间 约 7% ,主要 为 短期 急剧 的 干扰 事件 。 男 外 ,持续 时 间 为 4a 以 上 
四 县 不 超过 100hm?, 2000 年 之 前 持续 干扰 和 急剧 干扰 面积 相当 ,变化 比较 平缓 ;到 2000 年 之 
抽 者 都 呈现 波动 上 升 的 变化 趋势 。 在 选取 的 两 个 4km 的 样 方 中 ， 


TUR 


基于 影像 光谱 识别 以 及 通过 比 对 干扰 资料 的 可 视 化 验证 方法 表明 算法 结果 与 真实 地 表 的 解 译 信息 较 吻合 ,误差 约 为 0.1km 。 


F 


c 


定 基础 , 另 


一 方面 为 评价 森林 和 


关键 词 : 人 工 林 ;Landsat 时 间 序 列 ;LandTrendr; 森林 干扰 与 恢复 ;动态 监测 ;驱动 力 ;中 国 南 方 


书 长 时 间 序 列 遥 感 影像 进行 森林 干扰 的 自动 化 监测 十 分 必要 ,导出 的 定性 .定位 与 定量 信息 ,一 方面 为 可 持续 的 森林 经 营 黄 
产 力 与 森林 碳 储量 提供 有 效 的 数据 支撑 。 
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Abstract: Yearly Landsat imagery from 1986 to 2011 of a typical plantation region in Fogang County , Guangdong Province, 


southern China, was used as a case study. The pre-processing Landsat Ecosystem Disturbance and Adaptive Processing 


System ( LEDPAS) algorithm was implemented to generate standard surface reflectance images to construct a Landsat time 
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series stack ( LTSS). The LTSS was fed to the Landsat-based detection of trends in disturbance and recovery ( LandTrendr) 
algorithm to monitor long-term changes in plantation disturbance and recovery, followed by an intensive validation and a 
continuous 24 years change analyses on annual change, and disturbance amount and duration. Validations derived from two 
chosen sample plots of 4 km? indicated that the LandTrendr-based mapped disturbance results strongly agreed with those 
derived from the visual interpretation of the pre- and post-disturbance multispectral images and visualization of the local 
disturbance documents, with an error of 0.1 km?. Results indicated that the forest disturbances that occurred in Fogang 
County were relatively drastic. An annual disturbance of 1000 hm? was witnessed for most years of the study, and an annual 
disturbance of over 2000 hm? occurred in 1987, 2002, 2004, 2005, 2006, 2007, and 2009. Particularly, the disturbance 
of 1987 and 2007 exceeded 6000 hm”. In comparison to forest disturbance, forest recovery areas were relatively stable. 
Through a trend analysis of forest disturbance and recovery in Fogang County, forest disturbance and recovery areas mapped 
in the late 1980s through 1990s were less than those mapped after 2000, and the trend was lower than that after 2000. Since 
2000, the forest disturbance areas have gradually increased, with a slight increase in forest recovery, but the overall 
magnitudes of forest disturbance exceeded those of forest recovery. The area of forest disturbance with a duration of 1 year 
accounted for 3896, 2896 for a duration of 2 years, 2596 for a duration of 3 years, and 796 for a duration of 4 years; these 
disturbances were classified as abrupt and short-term disturbance events. Gradual forest disturbance and recovery events for 
a duration over 4 years existed, but the overall areas were less than 100 hm^/a, and were highly different from the areas of 
abrupt disturbance events. Prior to 2000, abrupt and gradual disturbance areas were almost equal, with a gentle change. 
After 2000, abrupt disturbance areas were greater than those of gradual disturbances, with a maximum of 2800 hm , and 
both abrupt and gradual disturbances showed an undulatory increasing trend. Based on the history and status of forest 
disturbances in Fogang County, the factors contributing to the environmental disturbance of forest plantations were analyzed 
to develop effective forest management strategies and countermeasures. The current study demonstrated the need to use dense 
time series images to map forest disturbance and recovery events in plantation forests. This approach could provide 
qualitative, locational, and quantitative forest. change results for the land use decision-makers and conservation 
communities, enabling the strategic development of sustainable forest management and provide effective data support to 


evaluate forest productivity and carbon storage. 


Key Words: plantation; landsat dense time series; LandTrendr; forest disturbance and recovery; dynamic monitoring; 


driving forces; Southern China 


干扰 是 森林 生态 系统 动态 变化 的 主要 驱动 力 ,干扰 的 历史 影响 林 分 的 生长 状态 ,不 同 干扰 的 类 别 .强度 和 
大 小 能 改变 林 分 物种 组 成 和 林 分 结构 …” 。 典 型 的 自然 干扰 ( 雨 雪 灾害 病虫害 等 ) 和 人 为 干扰 (采伐 、 土 地利 
用 变化 等 ) 以 及 干扰 后 更 新 事件 将 影响 森林 碳 汇 ” 。 不 同 景观 变化 下 ( 如 城市 化 气候 或 灾害 变化 及 森林 管 
理 活动 等 ) ,干扰 的 状态 受 不 同 程度 的 影响 。 和 森林 采伐 作为 人 为 干扰 及 景观 动态 的 重要 部 分 ,一 方面 影响 
着 森林 碳 储量 , 男 一 方面 它 与 森林 覆盖 度 、 物 种 分 布 地形 及 社会 经 济 因 素 等 密切 关联 '“ 。 自 从 1980 年 开始 ， 
大 尺度 的 造林 和 再 造林 计划 的 实施 使 得 中 国 具 有 了 最 大 的 人 工 林 分 布 面积 , 占 全 球 人 工 林 面积 的 1/4'”。 人 
工 林 的 存在 使 得 森林 碳 储 量 得 以 增加 中 。 在 南方 区 域 ,造林 和 再 造林 活动 是 区 域 碳 汇 变化 的 主要 原因 之 
一 。 人 工 速生 丰产 林 建 设 承 担 着 缓解 木材 需求 ,发 挥 改善 经 济 社会 .生态 效益 的 重要 作用 ,被 列 为 国家 产 
业 政 策 ,具有 举足轻重 的 作用 '"。 但 是 人 工 林 普 遍 存在 树种 单一 结构 和 功能 较 差 ,病虫害 频繁 发 生 的 特点 。 
另外 ,人 工 林 营造 时 树种 选择 主要 考虑 速生 丰产 \ 短 轮 伐 期 的 特性 "”" ,这 使 得 中 国 南方 区 域 的 森林 存在 快 
XE 广泛 的 时 空 变 化 。 发 展 科 学 可 靠 的 方法 来 监测 森林 干扰 或 者 恢复 时 空 模式 对 于 人 工 林 健 康 经 营 管理 及 其 
碳 核 算 意 义 深远 。 

然而 , 仅 凭 借 两 期 或 者 三 期 长 时 间 间 隔 的 影像 “中 往往 不 能 充分 捕捉 南方 人 工 林 频繁 发 生 的 采伐 及 伐 
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后 更 新 活动 。 近 年 来 基于 30m 分 辩 率 的 Landsat 卫星 遥感 技术 的 发 展 为 森林 干扰 的 监测 提供 了 数据 支撑 , 特 
别 是 时 间 序 列 的 遥感 数据 被 成 功用 于 森林 干扰 或 恢复 变化 监测 ”…: 。 但 是 如 何 能 快速 有 效 地 识别 出 不 同 持续 
时 间 和 地 域 分 布 的 森林 干扰 或 恢复 现象 仍然 是 项 技术 难题 。 尽 管 一 些 森 林 变 化 自动 化 分 析 算 法 :2 为 森林 
生态 系统 的 近 实 时 监测 提供 了 强 有 力 的 工具 ,但 是 很 少 有 方法 既 能 监测 短期 的 森林 变化 又 能 探测 持续 的 和 森林 
变化 现象 。 基 于 Landsat 时 间 序 列 数据 天 汪 并 充分 利用 不 同 干扰 类 型 的 时 间 和 强度 特性 识别 干扰 或 恢复 的 
技术 ,如 LandTrendr 算法 5 ,已 被 成 功用 于 西北 太平 洋 区 域 森林 干扰 的 监测 ”而 且 利 用 1986 年 起 至 
当前 ,近乎 连续 的 时 间 记 录 南 方 人 工 林 干扰 或 者 恢复 分 布 的 研究 也 未 见报 道 。 所 以 ,基于 监测 所 得 综合 性 和 
连续 性 的 人 工 林 干扰 和 恢复 数据 记录 对 于 南方 区 域 森 林 管 理 和 政策 制定 .生态 系统 保护 和 重建 ,生物 多 样 性 
保护 和 碳 评价 有 重大 价值 。 

本 研究 主要 目的 基于 不 同 来 源 的 Landsat 长 时 间 序 列 遥感 影像 数据 堆栈 p122r043 (条 带 号 ) ,利用 
LandTrendr 算法 分 析 广 东 省 佛 冈 县 的 森林 干扰 与 恢复 的 时 空 模式 ,包括 森林 干扰 年 份 森林 干扰 量 、 和 森林 干扰 
持续 时 间 等 ,并 借助 现 有 影像 以 及 历史 数据 与 实地 调查 记录 ,验证 森林 干扰 与 恢复 的 制图 精度 ,分 析 森 林 变 化 
的 驱动 力 及 识别 典型 干扰 类 型 时 空 特性 及 其 关联 要 素 。 


1 研究 区 域 与 数据 


佛 冈 县 (图 1) 位 于 广东 省 中 部 ,珠江 三 角 洲 边缘 ,清远 市 东南 角 ,境内 低 山 丘陵、 谷地 .平原 交错 ,大 多 在 
海拔 300m 以 下 。 全 县 山地 丘陵 区 大 部 分 是 人 工 林 , 主要 树种 为 校 树 、 马 尾 松 .湿地 松 、 杉 木 . 荷 木 . 歼 戎 LUE 
人 竹 、 荔 枝 和 龙眼 等 果树 ,目前 全 县 森林 以 速生 丰产 的 商品 林 和 经 济 林 占 主导 (图 1) 。 本 研究 选取 了 1986— 
2011 年 间 24 期 Landsat TMZETM+ 有 适 感 影像 (红色 方形 框 为 pb122r043 一 景 数据 边界 ) 为 研究 数据 ,该 一 景 数据 
仅 履 盖 佛 冈 县 大 部 分 区 域 (以 下 均 简称 佛 冈 县 ) 作 为 后 续 研 究 区 (图 1)。 
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图 1 研究 区 地 理 位 置 
Fig.1 Geographical location of the study area 


本 研究 选取 的 24 期 Landsat TM/ETM -3 E [5 2c 8 T E E tb Jot e] x Je EA] JU ER v 0000] Ej RES ( USGS/ 
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EROS) rft (EA P fij fk EROS) ,其 余 购 买 于 中 国 科 学 院 对 地 观测 与 数字 地 球 中 心 ( 中 科 院 地 面 站 ,以 下 简称 
BJGS ) 或 者 从 地 理 空 间 数 据 云 处 免费 获得 。 影 像 获 取 日 期 的 选择 必须 满足 植物 生长 旺盛 的 6 一 9 月 份 ,在 一 
些 低 纬度 地 区 可 以 顺延 至 5—10 月 份 ,以 减少 物候 对 于 影像 光谱 识别 的 影响 。 在 云 量 选 择 上 ,由 于 获取 完全 
无 云 影像 有 一 定 难 度 , 尽 可 能 保证 较 少 云 量 , 在 一 个 时 间 序 列 内 要 保证 一 定 比例 的 质量 较 好 的 数据 。 表 1 描 
述 了 研究 所 用 的 Landsat TMZETM+ 遥 感 数据 的 相关 信息 。 广 东 省 林业 生态 规划 设计 院 提供 了 佛 内 县 2005 
年 一 2011 年 森林 资源 数据 库 和 2005 年 一 2006 年 的 小 班 数据 用 于 识别 森林 类 型 和 不 同 地 表 类 型 (图 1) 。 借 
助 Google Earth 提供 的 高 清晰 图 像 了 解 现 地 情形 ,并 根据 现 地 调查 和 查阅 当地 历史 资料 及 向 有 关 人 士 咨 询 发 
现 该 区 域 存在 大 量 森林 采伐 , 曾 出 现 一 些 灾害 现象 。 


表 1 研究 用 遥感 影像 数据 介绍 


Table 1 Description of remotely sensed images used in the current analysis 


MEN LL LANE M PEAN 
Uses Satellite Sensor Source quality date Satellite Sensor Source quality 

1986-11-03 Landsat 5 TM BJGS 高 1987-06-15 Landsat 5 TM EROS 低 
1988-11-08 Landsat 5 TM BJGS 高 1989-06-20 Landsat 5 TM EROS 低 
1991-04-07 Landsat 5 TM BJGS 3% zu 1992-7-30 Landsat 5 TM BJGS 3% zu 
1993-10-05 Landsat 5 TM EROS 1% zi: 1994-11-09 Landsat 5 TM BJGS 高 
1995-10-11 Landsat 5 TM EROS 10% zziz 1996-07-10 Landsat 5 TM EROS 10% KHT 
1997-11-01 Landsat 5 TM BJGS 5% zu 1998-08-16 Landsat 5 TM BJGS 15% zx 
1999-11-15 Landsat 7 ETM+ EROS 1296 云 覆 盖 。” 2000-08-21 Landsat 5 TM EROS 1% KHH 
2001-00-17 Landsat 7 ETM+ EROS 5% 77 2002-10-14 Landsat 5 TM EROS 5% KHM 
2004-10-03 Landsat 5 TM EROS 高 2005-07-18 Landsat 5 TM EROS 396 zu 
2006-09-23 Landsat 5 TM BJGS 3% zu 2007-04-19 Landsat 5 TM BJGS 5% zu 
2008-07-26 Landsat 5 TM BJGS 3% 云 覆盖 2009-10-17 Landsat 5 T™ EROS 3% zu 
2010-03-26 Landsat 5 TM BJGS 3% zu 2011-08-20 Landsat 5 TM BJGS 5% zu 


2 研究 方法 


2. 预 处 理 与 光谱 指数 设 定 
2.1.1 KARIE 

将 以 上 两 种 来 源 数据 基于 LEDAPS 进行 标准 化 处 理 构 建 LTSS ,该 方法 于 沈 文 娟 等 .中 中 被 具体 描述 。 
2.4.0. zx BASE dE s 

为 了 进一步 保证 影像 数据 光谱 识别 的 有 效 性 ,基于 需要 修复 的 目标 影像 和 无 云 污染 的 参照 影像 ,通过 参 
照 影像 整合 创建 一 个 平均 的 参照 影像 (参照 影像 的 数量 至 少 3 一 5 幅 , 确保 影像 清晰 或 者 云 掩 膜 效 果 好 ) , 然 
后 将 目标 影像 与 参照 影像 进行 比较 ,通过 人 为 设 定 的 阔 值 进行 云 和 阴影 的 掩 膜 。 云 掩 膜 主 要 借助 波段 1 进行 
BUELE , 当 像 元 值 小 于 阔 值 在 影像 上 显示 为 黄色 , 即 为 需要 掩 膜 的 云 覆 盖 区 域 ;阴影 掩 膜 主 要 利用 波段 4 和 
波段 5 进行 闵 值 设 定 ,黄色 显示 区 域 表 明 需 要 掩 膜 的 区 域 ( 图 2) 。 
2.1.3 ”光谱 指数 选取 

在 Kennedy 等 ' 站 的 研究 中 , 归 一 化 燃烧 指数 (NBR) 对 于 捕捉 干扰 事件 具备 最 大 敏感 性 , 比 起 归 一 化 植被 
指数 ( NDV D 和 继 帽 变化 湿度 指数 (Wetness) 而 言 有 较 高 干扰 匹配 的 可 能 性 。 另 在 Cohen £ (2010) WA 
中 ,表明 其 对 于 干扰 具备 较 好 的 解释 能 力 。 本 研究 中 使 用 NBR' 作为 算法 监测 的 指数 ,该 指数 通过 近 红 外 
(NIR) 和 短波 红外 (SWIR ) 的 差 值 与 和 值 的 比值 用 于 监测 干扰 。 

NBR = (NIR-SWIR)/(NIR+SWIR) (1) 
该 比值 中 ,NIR 反映 健康 绿色 植被 ,SWIR 反映 岩石 和 裸 土 。 健 康 的 森林 有 高 的 NIR 值 和 低 的 SWIR fE, 
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并 具备 高 的 NBR (E. SM T URS ,SWIR 值 增加 ,NBR 
值 表 现 为 更 低 。 然 后 ,时 间 序 列 堆栈 中 的 影像 数据 被 转 
化 成 NBR 值 , 与 云 掩 膜 数 据 进行 像 元 与 像 元 之 间 的 匹 
配 。 并 将 这 些 时 间 序 列 源 数据 用 于 算法 分 割 。 
2.2 时间 序列 分 割 算法 与 滤波 
LandTrendr 算法 的 核心 为 时 间 序 列 分 制 。 分 割 算 
法 使 用 时 间 序 列 中 的 每 个 像 元 来 识别 变化 的 NBR 值 ， 
并 保留 时 间 序 列 中 NBR 值 的 最 大 变化 ,对 应 的 小 的 变 
化 进行 滤波 操作 , 目的 是 增加 信 噪 比 。 具 体 参 数 设置 如 
表 2 所 示 。 分 割 的 第 一 步 是 由 顶点 的 年 份 决定 ,用 于 定 
义 分 割 的 结束 点 ,第 二 步 是 基于 这 些 顶 点 ,使 用 点 与 点 
或 者 回归 线 的 灵活 混合 得 到 最 好 的 直线 轨迹 ,从 分 割 算 
法 得 到 的 值 有 :每 年 的 源 数据 (每 年 对 应 像 元 最 好 的 无 图 2 云 、 阴 影 掩 膜 阔 值 设置 
Z NBR 值 ) .顶点 对 应 的 年 份 .每 年 的 拟 合 NBR 数据 。 Fig.2 Threshold settings for cloud and shadow screening 
(分 割 中 每 个 点 的 NBR 值 被 用 于 描述 轨迹 ) ,这 些 数据 
被 用 于 后 续 的 制图 算法 。 整 个 过 程 识别 出 NBR 值 的 3 种 特性 , 即 增加 (恢复 ) .下 降 ( 干 扰 ) 和 保持 不 变 。 其 
中 基于 地 表 数 据 使 用 回归 模型 转换 NBR 值 来 评估 植被 覆盖 的 比例 生成 干扰 量 值 ”。 考 虑 到 该 回归 模型 是 
基于 西北 太平 洋 区 域 不 同 植被 覆盖 与 密度 构建 ,本 研究 区 植被 组 合 和 密度 与 其 的 不 一 致 性 ,在 算法 参数 设置 
时 进行 了 相应 的 调整 以 适应 目前 的 研究 区 的 植被 特征 ( 表 2)。 干扰 持续 时 间 值 为 1 一 25a( 取决 于 影像 堆栈 
的 年 份 数 ) ,由 下 降 的 NBR 值 对 应 的 年 份 数 决定 。 


表 2 LandTrendr 处 理 所 需 参数 


Table 2 Parameters used in LandTrendr processing 


"e dk "€ 25 
4] 3| Segmentation ”光谱 指数 NBR 滤波 Filter la 植被 覆盖 损失 阅 值 10 
轨迹 分 割 的 最 大 数量 6 20a 植被 覆盖 损失 阔 值 3 
拟 合 的 最 大 p 值 0.05 干扰 前 覆盖 阔 值 20 
KEE 0.25 TEE K LEE 5 
顶点 参数 0.9 制图 Mapping 最 小 制图 单位 11 像 元 / 约 Ihm? 
最 优 模型 比例 0.75 短期 干扰 /长 期 干扰 闭 值 10a 


2.3 ”和 森林 干扰 制图 

基于 年 份 干扰 量 和 干扰 持续 时 间 影 像 创建 干扰 图 (图 3) 。 最 小 制图 单位 用 于 确保 干扰 的 斑 块 至 少 有 11 
个 像 元 , 约 为 1 hm 的 大 小 ,此 浆 值 可 以 捕 换 到 大 多 数 的 森林 采伐 活动 ,也 避免 了 干扰 事件 制图 面积 太 小 而 使 
得 大 范围 内 验证 困难 的 问题 ( 表 2) 。 生 成 的 干扰 图 有 3 个 图 层 , 即 干扰 发 生年 份 . 相 对 干扰 量 .干扰 持续 时 
间 。 当 出 现 超 过 一 种 干扰 的 情形 ,根据 斑 块 面积 与 相对 干扰 量 的 乘积 之 和 进行 干扰 斑 块 的 分 类 ,来 平衡 不 同 
干扰 量 事 件 。 获 得 最 高 值 的 斑 块 即 为 主要 干扰 ,第 2 高 值 的 为 次 要 干扰 ,第 3 则 为 第 3 干扰 ,但 第 3 种 情况 较 
少 而 不 被 考虑 。 最 后 使 用 人 工 解 译 获取 的 森林 训练 区 样本 对 影像 进行 随机 森林 分 类 “i 。 掩 膜 中 识别 的 森林 
像 元 在 整个 时 间 序列 中 始终 保持 森林 的 状态 ,最 终 的 干扰 图 滤波 去 除了 非 森 林 像 元 中 的 任何 干扰 。 本 研究 获 
取 了 佛 内 县 研究 区 域 变化 最 大 的 干扰 模式 .持续 4a 之 内 变化 最 大 的 干扰 模式 持续 4a 之 内 变化 次 要 大 的 干 
扰 模 式 持续 超过 4a 的 干扰 模式 。 由 于 该 算法 能 监测 急剧 或 持续 干扰 (恢复 ) 事 件 ,将 发 生 于 4a 内 的 干扰 定 
义 为 急剧 干扰 或 恢复 ,超过 4a 则 定义 为 发 生 了 持续 干扰 或 恢复 。 另 外 整个 算法 处 理 的 是 整 幅 影 像 ,最 后 用 研 
究 区 域 边界 切割 干扰 结果 影像 。 
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干扰 创建 : 
据 干 扰 量 和 面积 分 类 1. 主 要 干扰 年 份 : | SF | 
干扰 年 份 
干扰 量 急剧 干扰 
持续 时 间 


空间 滤波 : 
据 年 份 组 合 版 块 次 要 干扰 
最 小 制图 单元 Z AEFI 
持续 干扰 


于 扰 前 覆盖 3. 第 三 干扰 TA 


图 3 LandTrendr 分 割 算法 检测 森林 干扰 和 恢复 事件 流程 图 


Fig.3 Work flow for detecting forest disturbance and recovery events using LandTrendr segmentation algorithm 


2.4 森林 恢复 制图 
为 了 获取 和 干扰 后 的 森林 恢复 ,本 研究 计算 了 植被 恢复 的 两 种 指标 (绝对 和 相对 指标 ) 。 绝 对 指标 与 干扰 
后 的 恢复 有 关 ,表示 干扰 发 生 后 5a ZE NBR 值 的 变化 。 计 算 公式 如 下 : 
ANBR,,, = NBR geas -NBR ginea o (2) 
式 中 ,NBRs 和 NBR onea wo 分 别 表示 干扰 发 生 后 拟 合 的 Sa 时 间 内 NBR 值 和 顶点 定义 恢复 分 割 开始 的 NBR 
值 (末端 为 干扰 的 分 制 )。ANBR ,ww 定义 为 时 间 序 列 开始 之 前 干扰 引起 的 森林 损失 。 干扰 后 恢复 值 RI, 称 
为 恢复 指标 ,计算 公式 为 : 


ANBR ww 
I ANBR, 93) 
式 中 ,分 母 在 干扰 分 割 时 随 着 NBR 值 变 化 ,并 先 于 恢复 发 生 。 
2.5 ”和 森林 干扰 产品 验证 

基于 LandTrendr 算法 得 到 佛 风 县 森林 变化 最 大 的 干扰 年 份 图 ,由 于 变化 最 大 的 干扰 年 份 图 记录 的 是 像 
元 发 生 多 次 干扰 ,只 记录 其 中 变化 最 大 的 一 次 干扰 , 即 记录 在 该 区 域 发 生 的 最 主要 的 干扰 ,也 为 典型 森林 变化 
区 域 。 典 型 区 域 一 般 分 布 着 主要 干扰 和 次 要 干扰 。 主 要 干扰 ,如 皆 伐 、 林 分 蔡 换 事件 ( 火 ) ,干扰 后 的 像 元 在 


影像 上 显示 不 再 为 森林 像 元 ;次 要 干扰 其 干扰 后 像 元 显示 的 亮度 大 , 比 起 干扰 前 像 元 绿色 程度 下 降 ,但 是 始终 
保持 为 森林 像 元 。 通 常情 况 下 主要 干扰 能 够 反映 区 域 干扰 事件 的 发 生 ,影响 次 要 干扰 的 主要 为 一 些 低 干 扰 量 
事件 ,一 般 发 生 在 主要 干扰 事件 之 后 ,与 主要 干扰 保持 紧密 的 关系 ,在 一 定 程度 上 ,次 要 干扰 事件 可 以 忽略 。 
再 结合 主要 干扰 和 次 要 干扰 图 的 分 布 ,通过 光谱 可 视 化 方法 对 于 主要 干扰 和 次 要 干扰 的 发 生 进行 目 视 识别 。 
由 于 森林 恢复 一 般 需 要 经 过 较 长 的 过 程 ,获取 当地 验证 资料 不 容易 ,本 研究 暂 不 考虑 森林 恢复 的 验证 ,主要 验 
证 佛 冈 县 变化 最 大 的 干扰 年 份 图 。 

本 研究 选取 佛 冈 县 2 个 典型 的 区 域 ,使 用 受 干扰 年 份 对 应 的 Landsat 影像 验证 干扰 制图 (图 4) ,使 用 2km 
x2km 的 样 方 ,将 其 和 加 到 时 间 序 列 堆栈 中 的 每 个 影像 上 ,并 在 每 个 影像 对 应 的 区 域 提取 干扰 区 ,将 目 视 解 译 
得 到 的 干扰 区 徐 加 到 主要 干扰 年 份 图 上 ,对 比 真实 解 译 的 干扰 面积 与 算法 分 割 所 得 的 面积 ,确认 正确 识别 的 
精度 。 参 与 验证 的 资料 主要 来 自 于 对 现 有 24 期 Landsat 影像 的 目 视 解 译 及 含 小 班 区 域 的 野外 实地 调查 资料 ， 
同时 结合 Google Earth 进行 辅助 判别 。 


3 结果 与 分 析 


3.1 验证 森林 干扰 产品 
基于 影像 光谱 目 视 识别 及 对 照 进 行 干扰 制图 精度 的 验证 ,很 容易 地 获取 到 干扰 发 生 的 典型 区 域 。 通 过 解 
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4 使 用 Landsat 影像 可 视 化 验证 森林 干扰 年 份 图 


Fig.4 Visual validation of forest disturbance using Landsat images 


译 典 型 区 域 1 每 年 的 Landsat 影像 ,可 知 图 4 典型 区 域 1 的 主要 干扰 发 生 在 2006 年 ,选取 2002 年 和 2005 年 
的 影像 进行 比 对 ,发 现 其 他 年 份 也 有 一 定量 的 干扰 ,但 都 不 是 主要 干扰 。 对 比 目 视 解 译 的 面积 和 干扰 年 份 图 
中 的 面积 , 目 视 解 译 的 干扰 面积 为 0.8768km” ,主要 干扰 年 份 图 中 自动 干扰 制图 面积 为 0.7659 km 。 通 过 解 译 
典型 区 域 2 每 年 的 Landsat 影像 ,可 知 图 4 典型 区 域 1 的 主要 干扰 发 生 在 2008 年 ,选取 2002 , 2005 、2006、 
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2007 年 的 影像 进行 比 对 ,发 现 各 年 份 都 有 一 定量 的 干扰 发 生 ,2008 年 的 干扰 为 主要 干扰 ,与 算法 所 得 的 主要 
干扰 的 定义 相符 合 。 对 比 目 视 解 译 面积 和 干扰 年 份 图 中 的 面积 , 目 视 解 译 的 面积 为 0.3165 km? ,主要 干扰 年 
份 图 中 的 自动 干扰 制图 面积 为 0.2511 km, 经 过 以 上 的 研究 , 目 视 解 译 的 结果 与 LandTrendr 算法 所 得 滤波 干 
扰 图 像 的 结果 有 一 定 的 误差 , 约 为 0.1km 。 此 外 ,通过 影像 光谱 的 识别 ,也 能 辨别 干扰 分 布 的 区 域 。 


3.2 ”和 森林 最 大 干扰 和 恢复 变化 的 模式 分 析 
3.2.1 森林 最 大 干扰 和 恢复 变化 时 间 分 布 特征 

图 5 为 森林 变化 最 大 的 干扰 年 份 图 ,最 大 的 干扰 即 
像 元 发 生 多 次 干扰 ,只 显示 从 分 割 开始 到 结束 最 大 的 一 
次 干扰 ,图 中 水 体 . 非 森 林 、 持 续 森 林 已 被 掩 膜 , 其 余 不 
同 颜色 表示 不 同 的 干扰 年 份 分 布 。 本 研究 基于 像 元 尺 
度 , 刻 画 了 不 同年 份 的 森林 干扰 面积 以 及 森林 恢复 面积 
(图 6)。 发 现 佛 内 县 的 森林 干扰 较为 剧烈 ,一 般 都 在 
1000 hm2? 左 右 。 而 1987 , 2002 , 2004 , 2005 , 2006 , 2007 
和 2009 年 的 干扰 面积 均 超过 2000 hm ,其 中 1987, 
2007 年 两 年 的 干扰 面积 达到 6000 hm 以 上 。 相 比 和 森林 
干扰 的 变化 , 佛 冈 县 的 森林 恢复 面积 随时 间 的 变化 相对 
平稳 ,以 1987 年 为 最 大 ,达到 了 19000 hm 。 通 过 对 佛 
办 县 森林 干扰 和 恢复 面积 的 趋势 分 析 , 发 现 20 世纪 80 
年 代 末 到 90 年 代 和 森林 干扰 和 恢复 的 面积 基本 少 于 
2000 年 以 后 的 变化 面积 ,变化 趋势 比 2000 年 以 后 的 平 
缓 ;从 2000 年 开始 ,森林 干扰 面积 逐渐 上 升 ,变化 趋势 
高 于 森林 的 恢复 ,但 森林 的 恢复 面积 仍 有 所 提升 。 

将 干扰 的 像 元 进行 邻近 斑 块 组 合 ,并 确保 干扰 斑 块 
的 最 小 面积 为 1lhm ,用 于 计算 不 同 干扰 持续 时 间 的 面 
积 百 分 比 。 图 7 显示 佛 冈 县 的 森林 干扰 面积 随 着 持续 
时 间 总 体 呈 下 降 趋 势 , 即 干扰 事件 的 发 生 时 间 保 持 较 短 
的 年 份 ,主要 发 生 在 1 一 5a 范围 内 ,其 中 1 一 2a 的 干扰 
持续 时 间 最 多 ,发 生 在 la 的 干扰 面积 所 占 比例 最 高 , 约 
占 26% ,表明 佛 由 县 的 森林 干扰 主要 为 急剧 的 干扰 
事件 。 
3.2.2 ”森林 最 大 干扰 变化 空间 分 布 特征 
佛 冈 县 的 森林 干扰 存在 较 强 的 空间 分 布 特征 ,图 8 
为 该 县 森林 干扰 面积 随 海 拔高 度 变 化 的 分 布 图 。 从 图 
中 可 以 看 出 随 着 海拔 的 升 高 ,研究 区 的 干扰 面积 百分比 
呈现 明显 的 下 降 趋势 ,其 中 ,海拔 分 布 于 125 一 179m 之 


干扰 年 份 

1987 m9 1993 mw 1998 m 2004 EE 2009 
EH 1988 E 1994 E 1999 NN 2005 Em 2010 
EH 1989 Œ 1995 Æ 2000 NH 2006 — WB 2011 
EE 1991 1-3 1996 = 200] EE 2007 

E 1992 mmu 1997 (32002 Nm 2008 


图 5 佛 冈 县 森林 最 大 干扰 年 份 图 


Fig.5 Fogang county forest greatest disturbance map 


20000 -a- 干扰 
18000 -会 恢复 


面积 Area/hn? 
E 
E 
E 
E 


图 6 佛 冈 县 森林 干扰 和 恢复 面积 随时 间 变 化 特征 
Fig.6 Temporal features of forest disturbance and recovery 


areas detected over Fogang county 


间 的 干扰 分 布 面积 最 大 。 总 的 来 说 ,超过 5096 的 森林 干扰 面积 发 生 在 海拔 小 于 179m 的 低 海拔 区 域 ,而 当 海 
RKT 787m 之 后 ,森林 干扰 面积 占 人 研究 区 森林 干扰 总 面积 的 比例 仅 为 0.96%。 


3.3 ”和 森林 急剧 干扰 模式 分 析 


佛 冈 县 主要 的 急剧 干扰 事件 (主要 干扰 ) 和 次 要 的 急剧 干扰 事件 (次 要 干扰 ) 均 得 到 识别 ( 表 3)。 人 研究 表 


明 20 世纪 80 年 代 末 到 90 年 代 末 森林 的 主要 干扰 面积 明显 小 于 2000 年 以 后 ,说 明 急 剧 干 扰 在 2000 年 以 后 
发 生 较 为 频繁 ,2005 .2007、2009 三 年 主要 干扰 面积 在 2000hm 以 上 ,2004、2006、2008 年 主要 干扰 面积 也 在 
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7 佛 闪 县 森林 干扰 面积 随 持续 时 间 变 化 分 布 图 8 佛 冈 县 森林 干扰 面积 随 海拔 变化 分 布 
Fig.7 Fogang county forest disturbance area changing with Fig.8 Fogang county forest disturbance area changing with 
duration altitude 
RI 森林 主要 和 次 要 干扰 的 面积 
Table 3 Area of forest major and minor disturbance 
类 型 
TUS 1987 1988 1989 1991 1992 1993 1994 1995 1996 1997 1998 1999 
Disturbance type 
NK 47.61 16.29 0.36 36.54 | 26.82 30.24 36.54 17.55 27.81 25.2 94.86 237.33 
Major disturbance 
Pd EHe 2 
次 要 抗 /bm 6.39 0.6 0.09 0.9 1.44 6.21 
Minor disturbance 
干扰 类 型 
2000 2001 2002 2004 2005 2006 2007 2008 2009 2010 2011 
Disturbance type 
二 EB 2 
主要 Turis 40.86 54 626.76 1011.96 2757.33 1634.31 2153.52 927.9 2051.73 521.46 46.53 
Major disturbance 
次 要 干扰 /hm” 


0.81 0.36 1.53 3.87 0.27 1.08 5.94 14.67 0.18 


Minor disturbance 


1000 hm 左右 ,到 2010 年 以 后 逐渐 下 降 , 其 中 2005 年 
主要 干扰 面积 达到 2700hm^, mH. 2000 年 以 后 主要 干 
扰 的 发 生 有 明显 的 波动 趋势 ,与 佛 冈 县 速生 丰产 林 、 人 
工 林 短 轮 伐 期 的 特性 相符 ,表现 出 每 隔 两 年 基本 出 现 的 
一 次 波动 ,这 种 情况 在 2010 年 之 前 特别 明显 。 图 9 为 
急剧 干扰 持续 时 间 分 布 ,主要 为 发 生 在 1 一 4a 内 的 干扰 
事件 ,每 年 的 干扰 面积 比例 各 不 相同 ,其 中 , 约 38% 的 


干扰 持续 la 的 时 间 , 约 28% 的 干扰 持续 2a E EO TT 
25% 的 干扰 持续 3a 时 间 , 约 7% 的 干扰 持续 和 时间, 说 | o Eu 
明 短期 的 干扰 事件 发 生 较 多 。 


34 ”森林 持续 干扰 和 恢复 模式 分 析 Elo 佛 冈 县 森林 急剧 干扰 持续 时 间 
由 于 1987 年 的 相关 数据 异常 , 移 除 该 年 的 数据 。 Fig.9 Fogang county abrupt disturbance changing with duration 

另外 持续 干扰 和 恢复 的 定义 而 言 ,2007 年 之 后 没有 统 

计数 据 。 研 究 发 现 持续 的 森林 干扰 和 恢复 在 佛 冈 县 有 一 定量 的 分 布 (图 10) ,但 总 的 面积 不 大 ,不 超过 100 

hm ,2000 年 之 前 的 森林 持续 恢复 面积 较 持续 森林 干扰 大 ,到 2000 年 之 后 持续 干扰 的 面积 呈现 波动 上 涨 的 趋 
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势 ,2007 年 的 时 候 达 到 持续 干扰 和 恢复 的 最 大 值 ,分 别 
为 380 hm? ,220 hm 。 此 外 ,从 2005 年 开始 ,森林 的 持 
续 恢 复 面 积 增加 ,表明 在 这 段 时 间 森 林 履 盖 保 持 生 长 与 
更 新 的 状态 ,而 在 1989 年 也 有 森林 增加 的 状态 。 另 外 ， 
通过 比较 佛 冈 县 森林 持续 干扰 和 急剧 干扰 的 面积 ,发 现 
2000 年 之 前 持续 干扰 和 急剧 干扰 面积 相当 ,变化 比较 
平缓 ,到 2000 年 之 后 ,急剧 干扰 的 面积 远大 于 持续 干 
扰 , 最 高 约 达 2800 hm ,但 两 者 的 变化 趋势 都 呈现 波动 
上 升 (图 11)。 


4 ”讨论 与 结论 


41 驱动 力 分 析 

佛 交 县 属于 人 工 林 分 布 典型 区 域 ,分 布 着 大 量 的 人 
工 林 、 速 生 丰 产 林 , 其 森林 干扰 剧烈 的 原因 主要 由 于 一 
些 急 剧 的 干扰 事件 引起 ,如 人 工 林 的 采伐 、 炼 山 还 有 一 
些 沿 道路 和 河流 的 城市 建设 .恶劣 气候 有 影响 等 ,与 已 有 
森林 资源 调查 和 小 班 资料 的 记载 较为 符合 。20 世纪 
80 一 90 年 代 广 东 省 的 林业 政策 使 得 大 片 造林 的 出 现 ， 
干扰 相对 较 少 ' 中 。 然 而 ,1995 年 佛 冈 县 实施 速生 丰产 
核 树 营 造 计划 , 轮 伐 期 为 6al21 Ti 2000 年 以 后 随 着 人 
工 林 .速生 丰产 林 达 到 成 熟 期 ,大 量 的 林木 被 砍伐 。 同 
时 ,为 了 进一步 的 造林 ,造林 前 炼 山 整地 行为 也 催生 大 
量 的 森林 干扰 事件 ?1 。 并 且 随 着 2004 年 经 济 开始 
增 速 ,城市 化 建设 不 断 发 展 ,靠近 道路 和 河流 的 森林 不 
断 受 影响 ,使 得 森林 面积 急剧 减少 ,干扰 的 面积 呈现 明 
tob gk ss IK UU. 2005 年 到 2007 年 分 布 着 较为 广泛 
的 火烧 迹地 ,说明 人 为 炼 山 在 这 些 年 份 密集 出 现 。2005 
年 至 2011 年 间 , 用 材 林 、 速 生 丰 产 林 、 薪 炭 林 采伐 加 剧 ， 
形成 了 众多 采伐 迹地 。 其 中 2008 年 的 冰冻 雨 雪灾 害 造 
成 了 严重 的 森林 损毁 ,随后 的 林业 快速 恢复 政策 有 助 于 


~ 


年 份 Year 


0 50 100 150 200 250 300 350 
面积 Area/hm? 
图 10 佛 冈 县 森林 持续 干扰 和 恢复 的 面积 对 比 
Fig.10 Comparison of gradual forest disturbance areas and 


gradual forest recovery areas occurred in Fogang county 
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图 11 佛 冈 县 急剧 干扰 和 持续 干扰 的 面积 对 比 
Fig.11 Comparison of abrupt and gradual disturbance areas in 


Fogang county 


PARRE EI o EET, ATAR .速生 丰产 林 比 较 适 宜 种 植 于 中 低 海 拔 区 域 ,一 方面 易于 种 植 管理 , 另 一 方面 
易于 运输 管理 ,海拔 越 高 ,森林 受 干 扰 可 能 性 越 小 。1987 年 和 1989 年 干扰 相关 数据 出 现 异 常 值 ,由 于 图 像 的 


全 云 量 覆 盖 ,数据 质量 较 差 ,不 能 正确 地 识别 地 表 地 物 ,对 于 干扰 和 恢复 的 错 分 现象 明显 。 


然而 ,在 驱动 力 分 析 中 ,并 没有 考虑 到 影响 人 工 林 和 森林 干扰 和 恢复 所 有 可 能 的 自然 人 为 驱动 力 因子 ,而 
且 没有 建立 它们 之 间 定 量 的 关系 , 仅 从 定性 角度 进行 了 分 析 。 


42 算法 与 技术 分 析 
4.2.1 ” 预 处 理 与 算法 分 析 


LandTrendr 算法 在 监测 干扰 年 份 .干扰 量 的 基础 之 上 ,能 够 监测 特定 年 份 内 发 生 的 干扰 ,并 能 够 识别 主要 
与 次 要 干扰 的 分 布 ,特别 是 能 够 监测 急剧 干扰 和 持续 干扰 的 情况 。LandTrendr 获取 到 的 是 根据 干扰 量 和 面积 
确定 的 最 大 干扰 发 生 的 情况 ,已 经 略 去 了 一 些 细微 干扰 的 发 生 。 此 外 ,LandTrendr 算法 在 像 元 尺度 基础 上 进 
行 不 同类 型 的 滤波 ,保证 了 干扰 斑 块 面积 的 可 识别 性 ,有 助 于 监测 一 些 典 型 干扰 事件 。 该 算法 需要 建立 在 
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Landsat 时 间 序 列 数据 堆栈 基础 之 上 ,影像 数据 配 准 \ 质 量 的 好 坏 以 及 数据 的 格式 都 将 会 影响 后 续 算法 的 精 
度 。 本 研究 中 原始 数据 有 不 同 的 数据 来 源 , 不 同 数据 来 源 数据 的 处 理 与 存档 方式 不 一 致 ,在 进行 预 处 理 之 前 
必须 进行 标准 化 处 理 操作 ,而 非 标 准 化 数据 的 标准 化 过 程 必 将 存在 误差 ; 云 掩 膜 .算法 参数 设置 时 人 为 阔 值 设 
定 使 得 结果 误差 存在 ;算法 要 求 的 影像 日 期 分 布 在 7、8 月 份 ,由 于 不 可 能 获取 到 如 此 确切 日 期 的 数据 并 且 要 
保证 好 的 数据 质量 ,月 份 上 的 偏差 产生 物候 的 影响 ;影像 1 到 2 年 的 间隔 使 得 不 能 很 好 捕 换 缺失 年 份 的 信息 ， 
影响 干扰 年 际 监测 结果 。 本 研究 中 含有 者 干 数据 质 量 不 好 的 年 份 , 致 使 结果 的 误差 存在 ,如 出 现 1987 年 异常 
高 的 森林 恢复 面积 值 。 基 于 已 有 研究 中 所 证 明 的 最 优 指数 基础 上 选取 了 单一 指数 用 于 算法 监测 ,不 同 指数 对 
于 算法 运行 好 坏 的 贡献 未 知 ,或 者 不 同 指数 对 于 不 同 地 域 的 监测 效果 可 能 存在 差异 性 ,可 以 尝试 多 种 指数 的 
比较 ,获取 最 优 指 数 的 监测 。 
4.04.23 ”验证 方法 分 析 

国外 的 验证 技术 较为 成 熟 "”” 引 。 如 基于 TimesyncL "4 工具 的 LandTrendr 算法 验证 能 够 自动 化 获取 解 
译 结果 ,并 与 算法 分 割 结 果 做 比 对 ,考虑 该 方法 在 国内 使 用 缺乏 历史 地 面 数 据 和 和 森林 资源 数据 以 及 高 清 分 辩 
率 的 影像 ,未 将 其 在 本 研究 中 推广 。 但 基于 本 研究 中 的 验证 方法 所 得 的 结果 , 发现 干 扰 监测 的 结果 与 目 视 解 
译 的 结果 相 比 呈 偏 小 的 趋势 ,主要 考虑 到 该 算法 对 高 密度 针 叶 林 覆盖 区 域 的 森林 干扰 监测 更 为 敏感 ,对 于 本 
研究 中 非 完 全 针 叶 覆 羡 区 域 存在 一 定 的 限制 ,并 导致 结果 的 错 判 和 误 判 。 综 合 现 有 分 析 结 果 , 尽 管 存在 识 
别 上 的 不 足 , 监 测 结 果 仍 具 备 可 参照 性 。 在 以 后 国内 推广 使 用 中 需 将 本 研究 的 方法 形成 具备 相应 标准 的 体 
系 ,归档 所 需要 的 森林 干扰 历史 资料 ,收集 丰富 的 地 面 数据 辅助 验证 ,使 得 方法 应 用 所 得 结果 更 为 可 靠 。 
4.2.3 算法 与 技术 应 用 分 析 

森林 干扰 的 4 个 特征 , 即 持续 时 间 .强度 .频率 和 面积 ”1; 。 本 研究 未 能 实现 森林 干扰 各 项 特征 的 量化 与 
分 析 。 由 于 连续 年 份 的 干扰 类 型 与 生物 量 历史 统计 资料 的 缺乏 , 仅 能 监测 出 急剧 和 持续 干扰 发 生 的 持续 时 间 
和 面积 ,对 于 具体 的 干扰 类 型 强度 及 其 发 生 原因 定位 不 足 '” |。 
43 结论 

本 研究 尝试 探索 南方 人 工 林 区 域 森 林 干 扰 与 恢复 的 影像 自动 化 监测 技术 。 尽 管 LandTrendr 算法 在 国内 
的 初次 使 用 存 有 一 些 技术 与 数据 选择 困难 ,但 是 通过 整合 与 改进 预 处 理 技术 能 够 成 功 的 进行 干扰 制图 ,并 基 
于 量化 的 森林 变化 进行 驱动 力 定性 分 析 ,不 仅 能 够 提供 木林 干扰 与 恢复 的 历史 ,也 组 成 了 未 来 评价 木林 碳 储 
量 和 气候 变化 下 生物 多 样 性 保护 的 关键 部 分 ,并 为 国内 关于 森林 干扰 与 恢复 长 时 间 变 化 的 自动 化 分 析 提 供 了 
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